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(M : ABSTRACT 

We observed the reflection nebula (RN) NGC 7023, with the SH module, and the 
long-slit SL and LL modules, of the InfraRed Spectrograph (IRS) on Spitzer. We 
also present InfraRed Array Camera (IRAC) and Multiband Imaging Photometer for 
^ . Spitzer (MIPS) images of NGC 7023 at 3.6, 4.5, 8.0, and 24 /xm. We observe the 

aromatic emission features (AEFs) at 6.2, 7.7, 8.6, 1 1.3, and 12.7 /im, plus a wealth of 



CsJ ■ weaker features. We find new unidentified interstellar emission features at 6.7, 10.1, 



O 



O 



15.8, 17.4, and 19.0 /xm. Possible identifications include aromatic hydrocarbons or 



tJ" . nanoparticles of unknown mineralogy. We see variations in relative feature strengths, 

^ . central wavelengths, and feature widths, in the AEFs and weaker emission features, 

^ . depending on both distance from the star and nebular position (SB vs. NW). 



Subject headings: dust, extinction — infrared: ISM — ISM: individual (NGC 7023) 
— ISM: lines and bands — ISM: molecules — reflection nebulae 



1. Introduction 

The interstellar AEFs at 6.2, 7.7, 8.6, 11.3, and 12.7 /im characterize the mid-infrared (mid- 
IR) emission of the diffuse interstellar medium (ISM) of our own and other star-forming galaxies. 
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Duley & Williams (1981) first suggested that the AEFs were due to aromatic hydrocarbons. Leger 
& Puget (1984) and Allamandola, Tielens, & Barker (1985) attributed the AEFs to to polycyclic 
aromatic hydrocarbon (PAH) molecules with ~50 carbon atoms. We present here Spitzer (Werner 
et al. 2004) images and spectroscopy of NGC 7023, a RN with AEFs, illuminated by the Herbig 
Be star HD 200775, at a distance of 430 pc (van den Ancker et al. 1997). 



2. Observations 

We imaged NGC 7023 with IRAC (Fazio et al. 2004) in all four channels (3.6, 4.5, 5.8 and 

8.0 ^m) in two epochs (2003 Oct 10 and 2003 Dec 18). The images were reduced with the Spitzer 
Science Center (SSC) IRAC reduction pipeline, and combined with the SSC Mosaicer. This pro- 
cessing included dark subtraction, flat fielding, mux-bleed correction, flux calibration, correction 
of focal plane geometrical distortions and cosmic ray rejection. Figures 1 and 2 show our final 
images at 3.6, 4.5 and 8.0 yum. 

We obtained images at 24, 70, and 160 /xm with MIPS (Rieke et al. 2004) on Spitzer, in the 
scan map mode at medium rate. We covered 10' x 30' at all three wavelengths. We reduced the 
MIPS images using the MIPS Instrument Team Data Analysis Tool (Gordon et al. 2004). The 
reduction of the 24 /xm image (Fig. 2) followed preflight expectations, with a 10% uncertainty on 
the final absolute calibration. The 70 and 160 /xm images will be presented in a later paper. 

We measured the Short-High (SH) IRS (Houck et al. 2004) spectrum (9.9 - 19.4 ^m; R = 
A/AA = 600) of NGC 7023 on 2003 Sept 19. We obtained Long-Low (LL) long-slit IRS spectra 
(14.0 - 32.9 jj,m; R = 60 - 130) at two positions on 2003 Oct 1. These observations occurred 
during the checkout phase of Spitzer, resulting in larger than normal uncertainties in the SH and 
LL absolute positions of ±5". We observed the Short-Low (SL) long-slit IRS spectrum (5.1 - 15.1 
fim; = 60 - 130) of NGC 7023 on 2003 Oct 25. IRS was sufficiently checked out by this point, 
resulting in a much smaller absolute positioning for the SL data of ±1". 

The data were reduced and spectra extracted using Cornell IRS Spectroscopy Modeling Anal- 
ysis and Reduction Tool (Higdon et al. 2004). The spectra were calibrated by applying the spectrum 
and model template of point source calibrator a Lac to the intermediate unflat-fielded product of 
the SSC pipeline. We subtracted blank sky spectra only from the SH spectrum, observed as Spitzer 
cooled, to remove thermal emission from the then warm (~45 K) baffles of the telescope. All 
spectra were extracted from fixed pixel ranges over the slit, except for SH where we extracted the 
spectrum from over the entire aperture. 

Figure 2 shows that the SL slit intersects the filaments seen to the NW of HD 200775. Our 
primary LL slit intersects both the NW filaments and HD 200775 (Fig. 2). The SL slit crosses our 
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secondary LL slit at Position A, shown on Figure 2. We obtained the SH spectrum at Position B, 
also marked on Figure 2. 



3. Results 

We present our spectra of NGC 7023 in Figures 3,4, and 5. Table 1 gives the central wave- 
lengths and nebular positions at which we discovered new interstellar, spectrally resolved, emission 
features. 

TABLE 1 
New Infrared Emission Features 



Ac Positions 

(l^m) Detected'* 

6.7 Sl-3 

10.1 B, Sl-3 

15.8 A, L6-7 

17.4 A, B, L4-6 

19.0 Ll-5 



(a) Positions where features are detected. The nebular offset from HD 200775 of each position is 
labeled on Figure 3 (Positions A, B), Figure 4 (Positions SI - S8) and Figure 5 (Positions LI - 
L7). 

Figure 3a shows the combined SL and LL spectrum (5 - 33 /xm) at Position A. The absolute 
flux calibrations of SLl and SL2 are uncertain by 25%, yet the relative agreement between them 
is much better than this factor. The SL and LL spectra were not expected to match and we had 
to divide the LL2 and LLl spectra by factors of ~2 - 3, after correction for beam size, to get 
all spectral segments to match. We emphasize that this composite spectrum of Position A shows 
qualitatively the spectral features observed in this nebular region, and should not be used for quan- 
titative analysis. The brightest spectral emission features seen in the spectrum are the AEFs at 6.2, 
7.7, 8.6, 11.3, and 12.7 /xm. Weaker emission features are also present, both spectrally resolved 
emission features and unresolved H2 emission lines. 

Figure 3b shows our SH spectrum (9.9 - 19.4 /xm) at Position B. Moutou et al. (2000) and 
Van Kerckhoven et al. (2000) have observed this wavelength region in NGC 7023 at a similar 
resolution but a different spatial position, using the SWS spectrometer on ISO, at a lower signal- 
to-noise. Position A is near the H2-emitting NW filament, while Position B is in a region with little 
H2 emission (Lemaire et al. 1996). These two positions sample different physical conditions, yet 
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show similar spectra. The higher spectral resolution of our SH spectrum (Fig. 3b) allows us to 
easily distinguish between unresolved emission lines, such as the - S(l) H2 line at 17.03 /um, 
and resolved emission features, such as the adjacent 17.4 yum feature. Weak emission features that 
are barely spectrally resolved in Figure 3 a are clearly resolved in Figure 3b. 

One of the great strengths of Spitzer is the multiplex advantage, and the ability to discern 
spectral and spatial differences, of the long-slit spectrographs SL and LL. In Figure 4, we present 
the results from our long-slit SL observations. The AEFs are the brightest spectral features at all 
positions. We detect, in SL spectra closest to the star (SI - S4), the split between the 7.6 and 7.8 
^m features that mainly comprise the "7.7" fim AEF (Fig. 4). 

The 1 1 .0 /im feature is clearly separated from the 1 1.3 fim AEF in Figure 3b. In Figure 4, we 
can observe the 1 1 .0 /im feature blended with the 11.3 /im AEF, in SL spectra closest to the star 
(SI - S3). This blend, and its spatial variations, causes the 1 1.3 /im AEF to appear to decrease in 
FWHM and shift to longer wavelengths with increasing distance, d*, from the star. 

Figure 5 illustrates our LL2 spectra (14.0 -21.1 /xm) along the long-slit of IRS-LL. Figure 5, 
together with Figure 3, demonstrate that long ward of the five bright AEFs, there continues to be 
complex spectral structure. Most prominent in Figure 5 is an emission feature at 16.4 fim, NW of 
the star (L4 - L7). Fig. 5 also illustrates emission features NW of the star, at 15.8 /im, 17.4 //m, 
and 17.8 /im. These features to the NW are confirmed by the SL+LL spectrum of Position A (Fig. 
3a) and the SH spectrum of Position B (Fig. 3b). Figure 5 also shows an emission feature both 
SE and NW of the star at 19.0 /im. Figures 3, 4, and 5 yield abundant evidence of pure rotational 
lines of H2, marked on Figures 4 and 5. Fuente et al. (1999, 2000) have previously studied pure 
rotational lines of H2 in NGC 7023. 

Figures 3, 4, and 5 all show a faint but non-zero continuum, observed at 5 - 20 /im, in addition 
to spectral features and lines. This continuum emission is also detected spectroscopically at 2 - 4 
/im in NGC 7023 (Sellgren, Werner, & Dinerstein 1983; Martini, Sellgren, & DePoy 1999), along 
with a strong 3.3 /im AEF and plateau of emission at 3.2 - 3.6 fim. These spectra demonstrate 
clearly that the emission in the IRAC 3.6 fim filter (3.2 - 3.9 fim) in NGC 7023 (Fig. 1) is due to a 
mix of this continuum emission with the 3.3 fim AEF and its accompanying 3.2-3.6 /im plateau 
emission. 

Sellgren et al. (1983) and Sellgren (1984) proposed that the 2-5 fim continuum in NGC 
7023 and other similar RN is due to non-equilibrium thermal emission from tiny grains (~1 nm), 
stochastically heated to high temperatures (~ 1000 K) for a brief time by single stellar photons. The 
IRAC 4.5 fim filter covers 4.0 - 5.0 /im, a region in which no significant PAH features have been 
identified. We believe, therefore, that the emission within the IRAC 4.5 fim filter is completely 
due to this tiny grain continuum emission. The IRAC 8.0 fim filter covers 6.5 - 9.3 fim. Figures 
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3a and 4 clearly demonstrate that the IRAC 8.0 fj,m filter is dominated by AEF emission at 7.7 and 
8.6 fj,m. 

The 20 - 33 //m spectrum of Position A (Fig. 3a) shows a strong rise to longer wavelengths, 
producing the emission in the 24 /xm MIPS images (Fig. 2). IRAS 25 //m observations of the RN 
23 Tau (Castelaz, Sellgren, & Werner 1987) agree with predictions (Draine & Anderson 1985) that 
emission in this wavelength region is due to stochastically heated tiny grains. We expect this to be 
true ofNGC7023 as well. 



4. Discussion 

4.1. Discoveries of new emission features 

We have discovered new ISM emission features at 6.7, 10.1, 15.8, 17.4, and 19.0 //m in NGC 
7023. These features are spectrally resolved, so are not atomic or molecular emission lines. All 
five new features are observed in multiple positions with SH, SL or LL, giving confidence that 
they are not instrumental artifacts. The 10.1 /xm feature seen in SL, and the 17.4 /xm feature seen 
in LL, are both clearly detected in SH. These ISM spectral emission features are among the first 
spectroscopic discoveries made by Spitzer. 

Currently, these features have no identification. Possible identifications could include PAH 
bands, such as C-H out-of-plane bending modes for 10. 1 //m (Hony et al. 2001), or C-C-C bending 
modes for 15.8, 17.4, and 19.0 /xm (Van Kerckhoven et al. 2000). The spectral structure at 6 - 9 
fj,m is quite complicated, with only tentative identifications for well-known features (Peeters et 
al. 2002), making it difficult to predict whether the 6.7 //m feature might fit into the PAH model. 
Alternate possibilities, particularly at the longer wavelengths, include nanoparticles (~ 1 nm) of 
a specific mineral composition. For instance, Molster et al. (2001) identify a mixture of PAH 
species between 3 and 12 yum, and various crystalline and amorphous silicates beyond 18 ^m, in 
the spectrum of the planetary nebula NGC 6302. Other mineral species, such as simple metal 
oxides, or other combinations of abundant refractory materials, also remain to be explored. 

4.2. Spatial variations in spectral features 

We observe marked changes in spectral features across NGC 7023. The nebula SE of the star 
is where the new 19.0 /im feature is brightest. To the NW of HD 200775, the 19.0 fxm feature fades 
and the 16.4 fj,m feature becomes the brightest feature in our LL2 spectra. 

Our long-slit observations show that the relative strengths of various pairs of emission features 
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vary with d*. We consider feature pairs observed in the same order and module of IRS and find 
that the ratios 7.8 /im/7.6 fj,m, 7.4 /xm/7.6 /im, 11.3 //m/11.0 //m, and 11.3 ymllJ jira increase 
with increasing d^. 

One of the most marked spectral changes is an apparent weakening of the 8.6 jim AEF with 
increasing rf*. Uchida et al. (2000) and Cesarsky et al. (2000) observe this same phenomenon in 
other RN, and attribute it to a broadening of the 7.7 ^m AEF with increasing d^. We confirm that 
the 7.7 /im AEF is broader, and also find that its central wavelength increases with increasing d^. 
This suggests that the wing of the 7.8 jim feature, as it grows in strength relative to the 7.6 jim 
feature, steadily overwhelms the 8.6 /im AEF until it is barely visible. The 6.2 /im AEF, like the 
7.7 //m AEF, appears broader at larger d^. The equivalent width of the 6.2 /xm AEF also markedly 
decreases with increasing rf*. 

Another striking spectral variation with d^ is the changes in the 11-14 /im region. Close to 
the star, distinct features at 11.0, 11.3, 12.0, and 12.7 jim can be distinguished. Far from the star, 
this spectral region can only be fit by a blend of the 11.3 jira AEF and a broad bump of emission, 
centered at 12.5 jira and having a FWHM of 2.0 jira. This broad 12.5 jira feature increases with 
increasing d^, and is not observed close to the star. The same 11-14 jira behavior with increasing 
has been observed in the RN Ced 201 (Cesarsky et al. 2000). 

4.3. Imaging 

In the inner regions of the nebula, we find strong similarities between the IRAC 4.5 /xm 
continuum image and the IRAC 8.0 jim AEF image. The ring of emission between HD 200775 
and the NW filaments is real, and is observed in ground-based images of the 3.3 jim AEF (An & 
Sellgren 2003). 

On the largest scales (~6', or ~0.8 pc), our 4.5, 8.0, and 24 //m images show an hour-glass 
shape, containing little IR emission, which is clearly outlined by a narrow rim of IR emission. The 
filaments within ~1' of HD 200775 define the narrow waist of the hourglass shape. This spatial 
structure has been observed before, in 1-0 ^^CO, 2-1 ^^CO, and 3-2 ^^CO (Gerin et al. 1998; 
Fuente et al. 1998), at 10 - 20" resolution. Watt et al. (1986) and Fuente et al. (1998) have argued 
that the hourglass shape is the fossil remnant of a bipolar outflow from the Herbig Be star HD 
200775, implying that the large-scale mid-IR emission outlines the walls of the cavity produced 
by the earlier outflow. 

Note added in manuscript. — After this paper was completed, M. Jura of UCLA kindly pointed 
out to us that the positions and widths of the new 17.4 and 19.0 /xm features agree rather well with 
those of the two lower frequency transitions of the Ceo molecule (Frum et al. 1991). It is not possile 
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to establish detection of interstellar Ceo on the basis of the data in the present paper, but we will 
explore this potential identification with additional observations and analysis. Note, however, that 
Moutou et al. (1999) have set limits on the abundance of C60 and Cqq in NGC 7023 based on 
nondetection of the two higher frequency vibrational transitions between 7 and 9 /im. 

We thank the people whose dedicated work made the Spitzer Space Telescope a reality (see 
Werner et al. 2004). This work is based on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology under 
NASA contract 1407. Support for this work was provided by NASA through Contract Number 
1257184 issued by JPL/Caltech. KS thanks the JPL Center for Long-Wavelength Astrophysics for 
support. 
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Fig. 1.— IRAC images of the central 6' x 6' of the RN NGC 7023, at 4.5 //m (top) and 8.0 //m 
(bottom). NGC 7023 spectra show that the IRAC 4.5 /im filter is continuum emission and the 
IRAC 8.0 ^m filter is primarily AEF emission. Both observations were made with the 12 s High 
Dynamic Range mode, consisting of a 1.2 s exposure for observing bright field stars, followed by 
a 10.4 s exposure for detecting the low surface brightness nebular emission. Even in this mode, the 
central star, HD 200775, is saturated in both images. North is up; east is to the left. 



Fig. 2.— MIPS image of the central ~12' x -10' of NGC 7023, at 24 fim (top), and IRAC image 
of the central ~3' x —3' region, at 3.6 /xm (bottom). We illustrate how the IRS slits overlay differ- 
ent parts of the central regions of our 3.6 //m IRAC image, which is our highest spatial resolution 
Spitzer image. We extracted and combined SL2 and SLl spectra at eight spatial locations. Spec- 
trum SI (42" W 8" N) is closest to the star and spectrum S8 (38" W 58" N) is farthest from the star. 
We extracted LL2 spectra at seven spatial positions. Offsets range from 25" E 19" S (spectrum LI) 
to 49" W 36" N (spectrum L7). Spectrum L3 is HD 200775. Position A (the intersection of SL 
and LL) and Position B (the slit location for SH) are marked. North is up; east is to the left. 



Fig. 3. — (a) Combined 5-33 /xm normalized spectrum (R = 60 - 130) of Position A in NGC 
7023 (top). We combined SL2 (5.1 - 7.6 fim; red), SLl (7.5 - 15.1 fim; green), LL2 (14.0 - 21.1 
^m; magenta), and LLl (20.9 - 32.9 /xm; blue) spectra. No scaling was done between SL2 and 
SLl. LL2 and LLl were multiplied by factors of order ~2-3, after scaling by the beam size, in 
order to match SLl. (b) NormaUzed SH spectrum (9.9 - 19.4 //m; R = 600) of Position B in NGC 
7023 (bottom), extracted from the entire entrance slit (4"7 x 11"3). The spectrum is a composite 
of eleven orders, with substantial overlap between each. Relative uncertainty in the flux calibration 
between orders is ~7%. Adjacent orders are shown in contrasting colors. 



Fig. 4.— Long-slit SL (R = 60 - 130) spectra of NGC 7023. We extracted and combined SL2 
(second order; 5.1 to 7.6 /im) and SLl (first order; 7.5 to 15.1 fim) spectra at eight spatial locations, 
with a 7"2 x 3"6 box. No scaling was done; SLl and SL2 agree well at overlapping wavelengths. 
The absolute flux uncertainty in SL is 25%. The y—axis is flux density (Jy). The panel for each 
of the spectra SI (bottom), S2, S3, S4, S5, S6, S7, and S8 (top) is labeled with the offset of the 
spectrum, from HD 200775, in arcsec, and the wavelengths of the H2 lines 0-0 S(5) 6.91 //m, 0-0 
S(3) 9.66 fj,m, and 0-0 S(2) 12.28 //m are marked (vertical lines). 
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Fig. 5.— Long-slit LL2 (14.0 - 21.1 fim; R = 60- 130) spectra of NGC 7023, extracted at seven 
spatial positions, each with a 15'.'3 x 10'.'6 box. The axis is flux density (Jy). The absolute flux 
uncertainty in LL2 is 25%. The panel for each of the spectra LI (bottom), L2, L3, L4, L5, L6, and 
L7 ( top) is labeled with the offset of the spectrum, from HD 200775, in arcsec, and the wavelength 
of the H2 line 0-0 S(l) 17.03 //m is marked (vertical line). Spectrum L3 is HD 200775. 
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